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ABSTRACT 

We present the discovery of eight new radio pulsars located in the Large Magellanic 
Cloud (LMC). Five of these pulsars were found from reprocessing the Parkes Multi- 
beam Survey of the Magellanic Clouds, while the remaining three were from an ongoing 
new survey at Parkes with a high resolution data acquisition system. It is possible that 
these pulsars were missed in the earlier processing due to radio frequency interference, 
visual judgment, or the large number of candidates that must be analysed. One of 
these new pulsars has a dispersion measure of 273 pc cm" 3 , almost twice the highest 
previously known value, making it possibly the most distant LMC pulsar. In addition, 
we present the null result of a radio pulse search of an X-ray point source located in 
SNR J0047.2-7308 in the Small Magellanic Cloud (SMC). Although no millisecond 
pulsars have been found, these discoveries have increased the known rotation powered 
pulsar population in the LMC by more than 50%. Using the current sample of LMC 
pulsars, we used a Bayesian analysis to constrain the number of potentially observable 
pulsars in the LMC to within a 95% credible interval of 570001 3 qqoq. The new survey 
at Parkes is ^20% complete and it is expected to yield at most six millisecond pulsars 
in the LMC and SMC. Although it is very sensitive to short period pulsars, this new 
survey provides only a marginal increase in sensitivity to long periods. The limiting 
luminosity for this survey is 125 mJy kpc 2 for the LMC which covers the upper 10% 
of all known radio pulsars. The luminosity function for normal pulsars in the LMC 
is consistent with their counterparts in the Galactic disk. The maximum 1400 MHz 
radio luminosity for LMC pulsars is ~ 1000 mJy kpc 2 . 

Key words: pulsars: general — stars: neutron — Magellanic Clouds 



1 INTRODUCTION 

Our nearest galactic neighbours, the Large and Small Mag- 
ellanic Clouds (LMC and SMC), are irregular galaxies with 
a different star formation history than the Milky Way. 
They therefore represent very interesting targets for pul- 
sar searc hes. Previous pu l sar surveys of the Magellani c 
Clouds bvlMcCulloch et all il983h.lMcConnell et all |l99ll) . 
ICrawford et all l^Qulf ), and ^Manchester et al.l i|2006l ) have 
discovered 13 radio pulsars in the LMC. In addition, two 
rotation-powered X -ray pulsars have been dis covered (see 
ISeward et all (|l984T ) and iMarshall et al.1 (|l998h ^ leading to 
a total of fifteen known rotation-powered LMC pulsars. Five 
radio pulsars were also discovered in the SMC in the sur- 



veys bvlMcConnell et all i|l99ll ), ICrawford et al.1 lJ200ll) . and 
IManchester et al.l (J2006I) . yielding a total of twenty rotation 
powered pulsars in the Magellanic Clouds. 

Pu lsar surveys with the Park es 20-cm multibeam re- 
ceiver l|Stavelev-Smith et all Il99q) have been carried out 
along the pla ne of the Galaxy ([Manchester et al.ll200ll) . at 
interm ediate ([E dwards et al 2001 ) and high Galactic lati- 
tudes (|jacobv et all 12003 ; Burga v et al. 20061 ). and in the 
Magellanic Clouds l|Manchester et al1|200 6). These surveys 
have been very successful. Improved computer processing ca- 
pabilities and search algorithms have inspired a reprocessing 
of these various Parkes Multibeam Surveys, most notably 
the Galactic plane surve y, and have led to a large number 
of additional discoveries (|Keith et al.ll2009l ; iMickaliger et al.1 
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120121 ; lEatough et all [20131 ; iKnispel et al.ll2013h . In this pa- 
per, we present the discovery of five new radio pul s ars in 
the LMC from a reanalysis of the [Manchester et al.l (| 20061 ) 
survey data and three from a high resolution survey of the 
LMC which we are conducting at Parkes for a total of eight 
new radio pulsars. 

In Section [2j we introduce the searches that were per- 
formed. In Section [3l we explain the data reduction stage, 
including the dispersion measure (DM) , period (P) , and ac- 
celeration search ranges. Our results are then presented in 
Section [4] We discuss the implications of our results in Sec- 
tion \S\ and finally, in Section [IJl we draw our conclusions. 



2 OVERVIEW OF THE SEARCHES 

Five separate searches using both new and archival data 
were performed in order to detect new radio pulsars in the 
SMC and LMC. Data fro m the Parkes Multibea m Survey 
of the Magellanic Clouds ([Manchester et al.ll2006l ) were re- 
processed and searched using three different sets of search 
parameters. The first pass focused on low dispersion mea- 
sures (up to 300 pc cm -3 ) and used a broad, coarse accel- 
eration search. The second trial extended the DM range to 
800 pc cm -3 and had a very fine acceleration search. The 
third search used the archival data to target beams that 
contained high- mass X-ray binaries (HMXBs) in the LMC. 
These survey beams were then reprocessed using a very in- 
tensive acceleration search. New data were processed from 
a deep search of an X-ray point source in the SMC and a 
new high resolution survey of the LMC using the Parkes 
Telescope. 



orbiting similarly high-mass, non-deg enerate com panions: 
PSRs J0045 -7319 (itsel f in the SMC. iKaspi et al.l |l994h) . 
J1638-4725 [|Lvnell2008l), J1740-305 2 (|Stairs et al.ll200ll ), 
and B1259-63 l|johnston et al.lll992l ). Of these, only PSR 
B1259— 63 is a member of a Be/X-ray binary system. Is it the 
case that accretion and/or eclipsing from the donor prevents 
radio emission from being detectable in many cases? Per- 
haps, but some HMXBs may be detectable as radio pulsars 
during the non-accretion phases of their orbits, if their orbits 
are sufficiently wide. PSR B1259— 63 is one of these wide- 
orbit, eccentric systems, and it has been extensively stud- 
ied as both a radio pulsar and transient X- ray source since 
its d iscovery more than twenty years ago (jJohnston et al.l 
Il992f ). particularly at or near its periastron passages, which 
occur every 3.4 years. This system continues to provide ex- 
citing results using multiw a yelength observations (see , e.g., 
Chernvakoya et all l|2009l ). iKhangulvan et alj (|2012l ), and 



Kong et al] (|2012h S. 



HMXBs from the Magellani c Clouds High-Mass X-Ray 
Binaries Catalog (HEASARC) (|Liu et alj|2005l) were iden- 
tified and checked for coincidence with the existing beams 
in the MC radio survey. These beams were then selected to 
run a directed search for radio emission using the Parkes 
Multibeam Survey data. Of the 128 HMXBs in the cata- 
log, 36 had positions that were either not in a coincident 
beam or found near the edge of a beam (therefore yielding 
low sensitivity), 16 were found in beams with excessive ra- 
dio frequency interference (RFI), which left 76 candidates 
to search. Since it is likely that a pulsar's signal would be 
highly accelerated, these beams were processed separately 
with a very fine acceleration search. Appendix A shows a 
list of all 76 HMXBs that were searched. 



2.1 Parkes Multibeam Survey 

Archival data from the Parkes Multibeam Survey of the 
LMC and SMC were reprocessed and searched by each team. 
The original observations took place between May 2000 and 
November 2001. Using a centre frequency of 1374 MHz and 
a bandwidth of 288 MHz split into 96 channels, data were 
sampled every 1 ms for a total o bserving time of 8400 s per 
pointing ({Manchester et al.ll2006l ). Due to increased comput- 
ing power and search algorithms, allowing us to use lower 
signal-to-noise ratios and higher acceleration ranges (see 
Section [3}, we were sensitive to lo wer luminosity and more 
highly accelerated pulsars than the lManchester et al] (|2006l ) 
survey. 



2.2 High- mass X-ray Binaries 

The Magellani c Clouds have 128 catalogued HMXBs 
(jLiu et al.ll2005l ). and this provided a large sample of sources 
to search for corresponding radio pulsar signals with a uni- 
form sensitivity. Multiwavelength detection of a HMXB as a 
radio pulsar provides insight into the evolution of these sys- 
tems and into the phys ics of the compan i on sta r (see, e.g., 
iThompson et al] (| 19941 ) and iGrove et all (| 19951 ) pertaining 
to PSR B1259-63). 

Despite extensive study of these HMXB systems (both 
in the Galaxy and in the Magellanic Clouds) at X-ray wave- 
lengths as both persistent and as transient/bursting sources, 
only four radio pulsars have been found to date to be 



2.3 An SMC Point Source 

We conducted a deep targeted radio search of SNR 
J0047.2— 7308 which was centr ed on the X-ray poi nt source 
identified within the SNR (see iDickel et all feoOll )). These 
observations were conducted to look for radio pulsations 
from a putative pulsar at that location. Two separate 
search observations of duration 5.3 and 5.6 hours were 
conducted with the Parkes 64-m radio telescope in Febru- 
ary 2008. These observations used the centre beam of the 
multibeam receiver with a centre frequency of 1390 MHz 
and 0.5 MHz filters giv ing a total bandwidth of 256 MHz 
({Manchester et alj|200ll ). The observing setup and the ob- 
servations themselves we re identical to the sear ch observa- 
tion reported in detail bv lCrawford et alj (J2009J) which tar- 
geted XTE J0103— 728, with the exception of the integration 
times and the sampling time used. In the radio observations 
reported here, a sampling time of 80 fis was employed to 
preserve sensitivity to millisecond pulsations from a young 
pulsar. 



2.4 High Resolution LMC Survey 

New surveys of the Large and Small Magellanic Clouds have 
recently been initiated using a new high resolution backend, 
the Berkeley-Parkes-Swinbur ne data recorder (BPSR), lo- 
cated at the Parkes Telescope (IKeith et alj|2010f ). These sur- 
veys were performed using the multibeam receiver at Parkes 
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with a centre frequency of 1374 MHz and a bandwidth of 
340 MHz. Each pointing was observed for 8600 s and the 
data were sampled every 64 /is (see Ridley & Lorimer, in 
preparation) . These surveys of the LMC and SMC began in 
May 2009 and we plan to continue them until 2016. They 
are ~ 20% complete with 49 out of 209 pointings observed 
and we report on the results to date. 



3 DATA REDUCTION 

The Parkes Multibeam data were searched multiple times 
with different search parameters. We first used the SIG- 
PROC processing packagqj to search 60 DMs from — 
300 pc cm -3 using a minimum signal-to-noise ratio (S/N) 
of 9.0, and this search was sensitive to pulsars with peri- 
ods from 10 ms to 10 s. A time-domain acceleration search, 
as implemented in the SIGPROC program seek, was per- 
formed over the range ±50 m/s 2 using intervals of 1 m/s 2 . 
Most pu lsars in binary sys tems have accelerations up to 
5 m/s 2 (|Camilo et al.ll2000l ). however, some, s uch as PSR 
B1744-24A with a=33 m/s 2 (|Lvne et al.lll990l ), can be sig- 
nificantly greater. To minimize the drifting of the signal 
across multiple Fourier bins due to the trial acceleration 
differing from the true acceleration, we chose a step size so 
that Aa = cP/T 2 bs . With the 8400 s observations of the 



Parkes Multibeam Survey, the acceleration search was fully 
sensitive to all pulsars in binary systems having spin periods 
greater than 235 ms. The data were also searched for single 
bursts of radio emission having a pulse width between 64 /j,s 
and 65 ms. A minimum flux density threshold for a single 
pulse search can be calculated (see ICordes fc McLaughlin! 
12003ft as 



Omin — 



(S/N), 



\v 



ir 



n p A/' 



(1) 



where W is pulse width, S sys is the system equivalent flux 
density (34 Jy averaged across all beams of the multibeam 
receiver), n p is the number of polarisations summed (2 for 
our case), and A/ is the receiver bandwidth. This made 
the search sensitive to giant bursts (of width 65 ms) having 
fluxes greater than 54 mjy. 

A deeper search was also performed using the SIGPROC 
processing package and searched 160 DMs over a range of 
- 800 pc cm -3 . The Fourier S/N threshold used was 5.0, 
and all candidates above that were folded. However, after 
using the prestoj program prepfold, only candidates with 
a resulting prepfold sigma 8.0 and greater were inspected. 
This allowed us to detect a candidate that might appear 
weak in the FFT but was then enhanced in the fold. Since 
sigproc conducts an incoherent sum of the harmonics in the 
Fourier search, and the folded data are coherently summed, 
a pulsar candidate can improve its signal-to-noise ratio sig- 
nificantly in the folded data and reveal an otherwise invisi- 
ble pulsar. This method is, however, more computationally- 
demanding. An acceleration search was also performed with 
a range between ±10 m/s 2 using intervals of 0.1 m/s 2 which 



1 http://sigproc.sourceforge.net 

2 http://www.cv.nrao.edu/ sransom/presto 



was fully sensitive to all pulsars with spin periods greater 
than 23.5 ms. 

Survey beams containing the 76 HMXBs were searched 
using periodicity, single pulse and acceleration searches. 
DMs up to 800 pc cm -3 (160 different DMs) were searched 
with candidates having a S/N greater than 5.0 being con- 
sidered. Accelerations between ±20 m/s 2 with step sizes of 
0.01 m/s 2 were tried, making this search fully sensitive to 
all binary pulsars having spin periods greater than 2.35 ms. 

The pointed search observations of SNR J0047.2-7308 
were analyzed to look for a radio periodicity from a possi- 
ble pulsar. After excision of RFI via the presto program 
rf if ind, each of the two data sets were dedispersed using 
a range of DMs from to 800 pc cm -3 , which easily en- 
comp assed the expected DM range for pulsars in the SMC 
IjManchester et al.ll2006l ) . Each separate dedispersed time se- 
ries was Fourier transformed, and the resulting spectra were 
filtered and harmonically summed, then searched for can- 
didate periodicities. Candidate signals were then checked 
by dedispersing and folding the data at DMs and periods 
near the candidate values. Each data set was also searched 
for dispersed impulsive signatures in case the source was 
an intermittent pulsar or rotating radio transient source 
IjMcLaughlin et al.l 120061) . The analysis closely followed the 
Fourier an d single pulse search a nalysis of similar data de- 
scribed bv lCrawford et al.l (|2009f l- No promising candidates 
were found in the Fourier or single pulse search of either 
data set. 



4 RESULTS 

Four pulsars were discovered in the reprocessing of the 
Parkes Multibeam Survey of the Magellanic Clouds data, 
while one pulsar was found in the search of HMXBs. An 
additional three pulsars were found in the initial process- 
ing of the high resolution LMC survey data. These pulsars 
were confirmed in separate observations at the Parkes Ra- 
dio Telescope between 14 Jan 2011 and 24 June 2012. Some 
general properties of these eight new radio pulsars can be 
found in Table [T] while their time-resolved and integrated 
pulse profiles are shown in Figure [l] 



4.1 Properties of 8 New Pulsars 

These 8 new LMC pulsars encompass a wide range of periods 
and dispersion measures. The periods range from 1.133 s for 
J0458-67 down to 112 ms for J0537-69. The 112 ms pulsar 
is the second shortest period known for a radio pulsar in 
the LMC after PSR B0540-69. It is also possible that it 
could be a young pulsar. Further timing observations are 
now planned for this and the other new discoveries. 

The range of dispersion measures of these pulsars is sig- 
nificant. Three of the newly discovered have DMs less than 
100 pc cm -3 . However, it is likely that they are all associ- 
ated with the LMC since their DMs are all >70 pc cm -3 , 
higher than the four lowest DM pulsars in the LMC, and 
well above the expected Gala ctic contribution to the DM in 
this direction (~25 pc cm~ 3 ; ICordes fc Laziol (120021 ) ). PSR 
J0537-69, on the other hand, has a DM of 273 pc cm" 3 . 
This is by far the highest DM of any known Magellanic 
Cloud pulsar (almost twice as high as the highest previously 
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Figure 1. Discovery observations and integrated pulse profiles of the eight new LMC pulsars. For each pulsar, the grey scale shows pulse 
intensity as a function of pulse phase and observation time. The sum of these data result in the integrated profile shown at the top of 
each plot. Two pulse phases are shown in each case for clarity. 
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known dispersion measure of 146 pc cm -3 for PSR B0540), 
and indicates that it is likely the most distant LMC pulsar 
yet discovered. 

4.2 SNR J0047.2-7308 

No radio pulses were observed from the point-like source 
in the supernova remnant J0047.2— 7308. The flux density 
limits on radio pulsations from the radio search depend on 
the period and duty cycle of any pulsations emitted by the 
source. For values that are typical for young pulsars (a 5% 
duty cycle and a spin period that is between a few tens and 
a few hundreds of milliseconds), the 1400 MHz flux den- 
sity upper limit to pulsed emission from the blind Fourier 
search is ~ 30 /xJy. For a distance of 60 kpc to the SMC 
i|Hilditch et al.ll2005r ) and assuming a 1 sr beaming fraction 
for any radio emission, this corresponds to 1400 MHz ra- 
dio luminosity of ~ 100 mjy kpc 2 . The limits on single- 
pulse emission are comparable to those rep orted for the ra- 
dio s earch of XTE J0103-728 in the SMC ([Crawford et all 
l2009h . 



5 DISCUSSION 

The known LMC pulsar population has now increased to 23 
pulsars with the addition of these 8 newly discovered pulsars. 
This is a 50% increase in the total number of LMC pulsars. 
In this section we discuss the DMs and flux densities of our 
new pulsars as compared to the previously known pulsars. 
Table [2] shows the properties of all currently known LMC 
pulsars. 
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Figure 2. The 1400 MHz flux densities of all LMC pulsars are 
shown as a function of their pulse periods. The 14 previously 
known radio pulsars (PSR J0537— 6910 is not included) are 
represented by filled circles while the 8 new pulsars and the lower 
limits of their fluxes are denoted by stars. The dashed line corre- 
sponds to the minimum flux threshold of the Manchester et al.l 
(2006) survey and the solid line is the minimum flux threshold 
of the new high resolution SMC and LMC surveys as ca l culate d 
following the procedure described in iManchester et al.l (2001). 
The increased sensitivities of the new survey are due to the 
increase in bandwidth, decrease in sampling time, and increase in 
number of frequency channels. The two pulsars below the survey 
thresholds (PSRs B0540-69 and J0535-66) were discovered 
using longer integration times. 



5.1 Flux Densities 

A compari son of the flux den s ities o f the previously discov- 
ered (from IManchester et al.l (|200q )) and newly discovered 
pulsars is found in Figure [2] We note that our improvements 
in computational power and search algorithms have enabled 
us to detect some pulsars with flux densities higher than the 
previous survey's minimum flux. The flux values for the eight 
new pulsars were obtained via the radiometer equation, 



(S/N) 



G y/n p i int A/ V P 



\Y 



W 



(2) 



where G is the telescope gain, n p is the number of po- 
larisations summed, T sys is the system temperature of the 
telescope, and W and P are the pulse width and period. 
For these calculations we use G=0.735 K/Jy for the centre 
beam, G=0.690 K/Jy for the inner ring and G =0.581 K/Jy 
for the outer ring (see IManchester et al.l ((2001J)), n p =2, and 
r aya =24 K. The pulse width, W, is assumed to be 0.05 P. 
Since the pulsars' positions are only known to within the 
beam radius, these quoted values are actually lower limits 
on their flux measurements. 



5.2 Dispersion Measures 

Most of the pulsars in both the Large and Small Magellanic 
Clouds have dispersion measures within the range of 65 — 



150 pc cm" 3 l|Manchester et alJl200r3 ). PSR J0537-69 was 
discovered to have a DM of 273 pc cm -3 which makes it 
the most highly dispersed pulsar located in the direction 
of the LMC and probably the most distant. The DMs of 
all LMC pulsars (with the exception of the X-ray pulsar, 
J0537-6910) are plotted in Figure El 



5.3 The population of LMC pulsars 

Using the current sample of pulsars in the LMC 
we applied the Bay es ian techniques developed by 
IChennamangalam et al.l (120130 to constrain the num- 
ber of potentially detectable LMC pulsars. Here, we take 
prior knowledge of the flux densities of the LMC pulsars 
and the mean and standard de viation (n and a) of the 
log-n ormal luminosity function (|Faucher-Giguere &: Kaspil 
J2006f ). the use of which is justified below in Section [5.41 to 
determine the population size (N). We follow the procedure 
and u sed the code detailed in IChennamangalam et al.l 
Q2013J) to carry out this analysis. Since the diffuse radio 
flux of the LMC is not uniquely connected to the pulsar 
population, our analysis did not make use of the tech- 
niques described in Section 2.3 of Chennamangalam et 
al., where the diffuse flux of the millisecond pulsars is 
used to constrain globular cluster pulsar populations. The 
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Table 1. List of newly discovered LMC pulsars, including the barycentred period, DM, epoch of discovery, integration time, and discovery 
survey. PM refers to the Parkes Multibeam Survey of the Magellanic Clouds reprocessing, BPSR is the High Resolution Survey, and 
HMXB is the targeted search of HMXBs. Parentheses represent the uncertainty in the last digit quoted. Right ascension is given in units 
of hours, minutes, and seconds, and the declination is given in units of degrees and arcminutes. The uncertainty in both RA and DEC 
is ±7 arcminutes. The uncertainty in DM was determined using prepf old plots, which is part of the PRESTO package. 







RA 


DEC 


P 


DM 


EPOCH 


T int 


Survey 


S/N 


PSR 




(J2000) 


(J2000) 


(s) 


(pc cm" 3 ) 


(MJD) 


00 






J0456- 


-69 


04 56 30 


-69 10 


0.117073051(15) 


103(1) 


52038 


8300 


PM 


29.2 


J0457- 


-69 


04 57 02 


-69 46 


0.231390388(73) 


91(1) 


55383 


8600 


BPSR 


11.3 


J0458- 


-67 


04 58 59 


-67 43 


1.1339000(18) 


97(2) 


51810 


8300 


PM 


10.8 


J0521- 


-68 


05 21 44 


-68 35 


0.43342070(30) 


136(4) 


51871 


8300 


PM 


19.3 


J0532- 


-69 


05 32 04 


-69 46 


0.57459786(70) 


124(1) 


55420 


7920 


BPSR 


9.6 


J0535- 


-66 


05 35 40 


-66 52 


0.210524357(30) 


75(1) 


51393 


14400 


HMXB 


7.1 


J0537- 


-69 


05 37 43 


-69 21 


0.112613212(20) 


273(1) 


55420 


7200 


BPSR 


8.5 


J0542- 


-68 


05 42 35 


-68 16 


0.42518900(72) 


114(5) 


51975 


3960 


PM 


18.6 



Table 2. Complete list of all 23 known rotation powered LMC pulsars. 







P 


DM 


Si 400 


Discovery 


PSR Name 


(ms) 


(pc cm" 3 ) 


(mjy) 


Paper 


J0449- 


-7031 


479 


65 


0.14 


Manchester et al. (20061 


J0451- 


-67 


215 


45 


0.05 


Manchester et al. (2006) 


J0455- 


-6951 


320 


91 


0.25 


McConnell et al. (1991) 


J0456- 


-69 


117 


103 


>0.15 


This work 


J0456- 


-7031 


800 


100 


0.05 


Manchester et al. (2006) 


J0457- 


-69 


231 


91 


>0.05 


This work 


J0458- 


-67 


1133 


97 


>0.07 


This work 


J0502- 


-6617 


691 


68 


0.25 


McConnell et al. (1991) 


J0519- 


-6932 


263 


119 


0.32 


Manchester et al. (2006) 


J0521- 


-68 


433 


136 


>0.12 


This work 


J0522- 


-6847 
-6652 


671 
975 


126 
103 


0.19 
0.30 


Manchester et al. (2006) 


J0529- 


McCulloch et al. (1983) 


J0532- 


-6639 


612 


69 


0.08 


Manchester et al. (2006) 


J0532- 


-69 


574 


121 


>0.05 


This work 


J0534- 


-6703 
-66 


1817 
210 


94 

75 


0.08 
>0.03 


Manchester et al. (2006) 


J0535- 


This work 


J0535- 


-6935 


200 


89 


0.05 


Crawford et al. (2001) 


J0537- 


-69 
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Figure 3. The DMs of all LMC pulsars (with the exception of 
J0537— 6910, which has not been radio detected) are displayed 
as a function of spin period. Filled circles represent previously 
known pulsars and stars correspond to the new discoveries. 
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Figure 4. The posterior probability density function for the to- 
tal number of potentially observable pulsars. The hatched areas 
represent regions outside of the 95% credible interval. 



X-ray pulsar, J0537— 6910, does not have a measured radio 
flux and hence is not included in our analysis. Likewise, 
two pulsars (J0535-6935 and B0540-6919) were detected 
in targeted searches using much longer integration times 
than the surveys discussed here. They were also excluded 
yielding a total of 20 pulsars for this analysis. 

We assumed a dist ance of 49.97 ± 1.3 kpc to the LMC 
IjPietrzvriski et al.ll2013 j ) and a minimum survey sensitivity 
of 0.05 mjy. Foll owing Bovles et al.l l|201ll ). we used the re- 
sults of lRidlev fc Lorimeil ( 2010a]) to constrain the luminos- 
ity function priors to be flat in the ranges —1.19 < /i < 
— 1.04 and 0.91 < a < 0.98. This results in a posterior prob- 
ability density function (see Figure^]) for the total number of 
potentially observable pulsars with a median of 57000^3qqqq 
where the upper and lower limits are given by a 95% credible 
interval. 

In an effort to further explore the luminosity function 
and underlying pulsar population we also tried a wide range 
of priors for fi and a. In this analysis we allowed /i to vary 
between —2.0 and 0.5 and a to vary between 0.2 and 1.4 ac- 
cording to lBagchi et al.l (|2011f ). The resulting constraints on 
the number of pulsars and the luminosity function are sum- 
marized in Table [3] The constraints on the total number of 
pulsars for these wide priors are consistent with the narrow 
prior results, albeit with large uncertainties. The currently 
observed population of LMC pulsars is not large enough to 
uniquely constrain the shape of the pulsar luminosity func- 
tion. More sensitive surveys with future instruments would 
undoubtedly help to sample the fainter regions of the lumi- 
nosity function and result in tighter constraints. 



Table 3. Median values for N, fi, and a for our Bayesian analysis 
of the LMC pulsar population and luminosity function. For the 
narrow priors case the quoted values for fi and <r only represent 
the range of priors that were used in the analysis. 



WIDE 
NARROW 



N 

21000+ 71000 
ziuuu_ 20000 

57000+ 70000 
°' UUU -30000 



1.55+i;g 1.14 



+0.22 
' -0.28 



-1 12+ ' 07 Q5+ - 03 
liX -0.07 u - ad _0.04 



5.4 Luminosity Distribution 

The above analysis assumes the log-normal luminos- 
ity dis tributio n found from s tudies of Galactic pulsars 
(Faucher-Gigue re fc Kaspill2006h . An interesting question to 
address with the growing sample of pulsars in the LMC is 
whether their luminosity distribution is indeed consistent 
with the Galactic population. To quantify this, we approxi- 
mate the high- luminosity tail of the distribution as a power 
law, where the number of pulsars (N) as a function of lumi- 
nosity (L) can be written as follows: 



log N = a log L + C. 



(3) 



Here C is a constant and a r epresents the slope of the distri- 
bution IjLorimer et alJuOOq ). We obtain a slope of a = —1.2 
for Galactic pulsars with L > 30 mjy kpc , and a slope of 
a — —3.6 for LMC pulsars with L > 125 mjy kpc 2 (see 
Figure [SJ|. At first glance, then, it appears that the slope 
and maximum luminosity of the two samples are markedly 
different. 

To investigate this apparent discrepancy between the 
Galactic and LMC luminosity distributions we simulated 
the luminosities of pulsars found in the LMC and com- 
pared them with the observed distribution. We first gen- 



© 0000 RAS, MNRAS 000, 000-000 



8 Ridley et al. 



erated 20 pulsars using a power law model with a = — 1.2_| 
and plotted the resulting cumulative distribution function 
(CDF) along with the actual observed luminosities. We 
then simulated another set of 20 pulsars using the log- 
normal luminosity function valu es of fi — —1.1 and a = 0.9 
IjFaucher-Giguere fc Kaspil [2000) . We ran 4 trials of each 
simulation type and plotted the results in Figure [5] As can 
be seen, the luminosities of the simulated pulsars do not 
vary significantly from the observed luminosities. We con- 
clude that, given the present statistics, the luminosity func- 
tion of normal pulsars in the LMC is consistent with that of 
their counterparts in the Galaxy. Our assumption of the log- 
normal luminosity function in the previous section therefore 
appears to be completely justified by the present sample. We 
attribute the steeper slope of the observed LMC pulsars as 
a bias due to the smaller sample size. The higher maximum 
luminosity seen for Galactic pulsars may be due in part to 
significant distance uncertainties associated with the Cordes 
& Lazio electron density model which are not present in the 
LMC sample due t o our adoption of a common distance from 
IPietrzvriski et al] IJ20131 ). A cap on the luminosity distribu- 
tion, based on the LMC sample, seems to be approximately 
1000 mJy kpc 2 . 



5.5 Small Magellanic Cloud 

Our reprocessing also searched the areas of the sky popu- 
lated by the Small Magellanic Cloud. However, no new pul- 
sars were detected. The smaller size of the cloud inherently 
leads to a fewer number of stars being pres ent, hence the like- 
lihood of detecting a pulsar is diminished. iRidlev fc Lorimerl 
i|2010br ) estimated that there are 40% fewer pulsars in the 
SMC than in the LMC, or N S mc = 0.6 Nlmc- Further de- 
creasing our odds of a new discovery is the increased distance 
to the SMC of 60 kpc (JHilditch et alj2005h . Since S = Lxd 2 , 
the observed flux of an SMC pulsar as compared to an LMC 
pulsar is given by 5*smc = Slmcm^Ip) , meaning that the 
observed flux is reduced by 35%. The five currently known 
SMC pulsars are the most distant pulsars confirmed to date. 



6 CONCLUSIONS 

With the 8 new pulsars presented in this work there are now 
23 known spin-powered pulsars known in the LMC, repre- 
senting a 50% increase in the population. Based on these 
results, using a Bayesian analysis, we constrain the total 
population of potentially observable pulsars in the LMC to 
be 57000t3°ooo ( 95% credible interval). The luminosity func- 
tion of LMC pulsars is consistent with that of normal pulsars 
in the Galaxy. Our results suggest that the maximum 1400- 
MHz radio luminosity is 1000 mJy kpc 2 . Further work is be- 
ing done to complete the high resolution BPSR survey of the 
LMC with the goal of establishing a deep and comprehensive 
survey of the entire galaxy. We hope to further increase the 
total number of known LMC pulsars an d find the first ex- 
tragala ctic millisecond pulsars (see, e.g., IRidlev fc Lorimerl 
l|2010br i). 



3 Similar re sults are found if we a dopt a flatter slope of a = —0.7 
as found bv lLorimer et al] (2006). 



We are also conducting a BPSR survey of the SMC 
(Ridley & Lorimer, in preparation) to find fainter and more 
rapidly rotating pulsars in that galaxy. After observing 12 
out of 73 planned pointings, no new pulsars have been de- 
tected, but a similar goal exists of performing an exhaustive 
search for pulsars in the SMC. 
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Table Al. Catalog of 76 HMXBs in the LMC and SMC searched for radio pulsations and dispersed radio bursts. The X-r ay type refers 
to whether the source is a pulsar (P), transient X-ray source (T), or ultra-soft X-ray spectrum (U) (see iLiu et al.l I I2005I) ). The RA is 
given in units of hours, minutes, and seconds, and the DEC is given in units of degrees, arcminutes, and arcseconds. 



RA (.12000) DEC (J2000) Orbital Period (d) X-ray Type Spect Type 



Spin Period (a) 



RX .10032.9-7348 

RX J0041. 2-7306 

AX J0042. 0-7344 

RX J0045.6-7313 

RX J0047. 3-7312 

AX J0048. 2-7309 

RX J0048. 5-7302 

RX J0049. 2-7311 

RX J0049.5-7310 

RX J0049. 7-7323 

XTE J0050-732#l 

XTE J0050-732#2 

RX J0050. 7-7316 

RX J0050. 9-7310 

XTE SMC46 

RX J0051. 3-7250 

AX J0051. 4-7227 

AX J0051.6-7302 

XTE J0051-727 

RX J0051.9-7255 

XTE J0052-725 

2E0051. 1-7304 

RX J0052. 9-7158 

CXOU J005323. 8-722715 

XTE SMC95 
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H 0053-739 
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